T - -  ■■ 

AD-A159  252  DETERMINATION  OF  NECESSARV  COMMANDS  FOR  THE  NOTION  OF  A  1/1 
ROCKET  ON  A  GIVEN.  .  (U)  FOREIGN  TECHNOLOGV  DIV 
HRIGHT-PATTERSON  AFB  OH  A  CODOBAN  27  AUG  85 
UNCLASSIFIED  FTD-ID<RS)T-1271-84  F/G  19/7  NL 


Wff  FILE 


FTD-ID(RS)T-1271-84 


FOREIGN  TECHNOLOGY  DIVISION 


DETERMINATION  OF  NECESSARY  COMMANDS  FOR  THE  MOTION  OF  A  ROCKET  ON  A 
GIVEN  LINEAR  TRAJECTORY  IN  THE  ATMOSPHERE 

by 

A.  Codoban 


dt \c 

elect fc 


Approved  for  public  release 
distribution  unlimited. 


EDITED  TRANSLATION 


FTD-ID(RS)T-1271-84  27  August  1985 

MICROFICHE  NR:  FTD-8 5-C-0007  3  4 

DETERMINATION  OF  NECESSARY  COMMANDS  FOR  THE  MOTION  OF 
A  ROCKET  ON  A  GIVEN  LINEAR  TRAJECTORY  IN  THE  ATMOSPHERE 

By:  A.  Codoban 

English  pages:  20 

Source:  Studii  si  Cercetari  de  Mecanica  Aplicata, 

Vol .  42,  Nr.  3,  May-June  1983,  pp.  223-239 

Country  of  origin:  Romania 
Translated  by:  SCITRAN 

F33657-81-D-0263 
Re  que  s  t  e  r :  FTD /TQTA 

Approved  for  public  release;  distribution  unlimited. 


THIS  TRANSLATION  IS  A  RENDITION  OF  THE  ORIGI¬ 
NAL  FOREIGN  TEXT  WITHOUT  ANY  ANALYTICAL  OR 
EDITORIAL  COMMENT.  STATEMENTS  OR  THEORIES 

PREPARED  BY: 

ADVOCATEDOR  IMPLIED  ARE  THOSE  OF  THE  SOURCE 
ANOOO  NOT  NECESSARILY  REFLECT  THE  POSITION 

translation  division 

OR  OPINION  OF  THE  FOREIGN  TECHNOLOGY  DU 

FOREIGN  TECHNOLOGY  DIVISION 

VISION. 

WP-AFB,  OHIO. 

e  -o 


DETERMINATION  OF  NECESSARY  COMMANDS  FOR  THE  MOTION 

OF  A  ROCKET  ON  A  GIVEN  LINEAR  TRAJECTORY  IN  THE  ATMOSPHERE 


By  Alexandru  Codoban 
Polytechnic  Institute  of  Bucharest 


The  variation  in  time  of  the  attack  angle  of  the  control  surfaces 
of  a  rocket  so  that  it  is  moving  on  a  given  linear  trajectory  in 
the  atmosphere  is  determined  in  this  paper.  Variation  of  aero¬ 
dynamic  coefficients  with  the  Mach  number  is  considered  as  well 
as  the  variation  in  time  (in  the  fuel  combustion  stage)  of  the 
moments  of  inertia  and  the  position  of  the  center  of  mass  of  the 
rocket.  Air  density  and  sound  velocity  are  assumed  variable  with 
the  altitude.  A  numerical  calculation  example  solved  with  the 
algorithm  presented  in  this  paper  is  shown  in  the  end. 


1.  MAIN  NOTATIONS  USED 


Oxy  is  the  reference  system  stationary  with  respect  to  the  Earth, 
with  the  Ox  axis  along  with  the  linear  trajectory  required  (fig.  1) 

-  the  angle  of  the  linear  trajectory  in  the  horizontal 

0  direction 

-  the  time  measured  from  the  moment  the  fuel  combustion 

t  begins 


v  —  -5l£.  —  i  ~  velocity  of  the  center  of  mass  of  the  rocket 

d t 

-  acceleration  of  gravity 

-  the  flight  altitude  of  the  rocket 

-  the  density  of  the  atmosphere  at  flight  altitude 

-  sound  velocity  in  atmosphere  at  flight  altitude 

M  -  Mach  number  for  the  flight 


/ 


O'l  T) 


-  reference  system  stationary  with  respect  to  the 
body  of  the  rocket  and  whose  axis  0'£ coincides 
with  the  axis  of  symmetry  of  the  rocket  (Fig.  1)_, 

-  mass  of  the  rocket,  variable  during  the  fuel  com¬ 
bustion  stage 

-  moment  of  inertia  of  the  rocket  with  respect  to 
the  rolling  axis  (which  passes  through  the  center 
of  mass) ,  variable  while  the  engine  is  running 

-  flow  of  gases  in  the  engine,  adimensional 

-  abscissa  of  the  center  of  mass  of  the  rocket  in  the 
system  O'  i;  r(; 

-  abscissa,  in  the  system  of  the  point  F  with 

respect  to  which  the  momentary  aerodynamic  coeffi¬ 
cients  are  calculated 
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-  traction  force  of  the  engine  of  the  rocket 

-  the  time  the  fuel  is  burnt  (including  t.  ) 

-  time  of  transition  from  active  stage  witn  engine 
running  to  passive  stage 

-  reference  surface  to  which  the  aerodynamic  coeffi¬ 
cients  are  referred 

-  the  reference  length  to  which  the  coefficient  of 
the  moment  of  rolling  is  referred 

-  the  angle  of  incidence  of  the  rocket  (the  angle 
between  the  direction  of  the  velocity  and  the  axis 
of  symmetry  of  the  body  of  the  rocket)  ;Vig .  1); 

-  the  angle  of  attack  of  the  surfaces  of  command  (fig.  1) 

-  angular  rolling  velocity  of  the  rocket 

-  angular  rolling  acceleration  of  the  rocket 

-  rolling  velocity  of  the  rocket,  adimensional 

-  aerodynamic  lift  force  of  the  rocket 

-  aerodynamic  drag  force  of  the  rocket 

-  aerodynamic  rolling  moment  of  the  rocket  with  respect 
to  point  F 

-  aerodynamic  lift  coefficient 

-  aerodynamic  drag  coefficient 

-  aerodynamic  drag  coefficient  at  zero  incidence 

-  base  pressure  drag  coefficient  of  the  rocket 

-  constant  in  the  formula  of  CD 

-  rolling  moment  coefficient  with  respect  to  point  F 
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derivatives  of  aerodynamic  lift  and  rolling 
moment  coefficients  with  respect  to  the  in¬ 
cidence  of  the  rocket 


a  »  d  =d0» 
dp  dp 

ci  =  ^ 

d$ 


derivatives  of  aerodynamic  lift  and  rolling 
moment  coefficients  with  respect  to  the  angle 
of  attack  of  the  command  surfaces 
derivative  of  the  rolling  moment  coefficient 
with  respect  to  the  adimensional  rolling  velocity 
(coefficient  of  rolling  reduction) 


DERIVATION  SIGNS  AND  INDEXES 


(* 


( 

( 


derivative  with  respect  to  time 

value  at  initial  moment  t  =  o 

o 

value  at  final  moment  t^ 


2.  INTRODUCTION 


Solving  a  problem  of  the  dynamics  of  the  flight  of  a  rocket  in 
the  atmosphere  is  the  scope  of  this  paper,  a  problem  often  met  in 
certain  practical  cases  like  the  determination  of  the  variation  law 
of  the  attack  angle  of  the  command  surfaces  so  that  the  rocket  will 
follow  a  given  linear  trajectory  in  the  atmosphere  which  it  must 
pursue  both  in  the  active  stage  of  the  flight  (with  the  engine 
running)  and  passive  stage  (without  traction  force) .  The 
variation  of  the  aerodynamic  coefficients  with  the  Mach  number  is 
being  evaluated  as  well  as  the  variation  in  time,  with  respect  to 
fuel  consumption,  of  the  inertial  moments  and  of  the  position  of 
the  center  of  mass  of  the  rocket.  These  are  considered  as 
numerical  functions  with  respect  to  the  Mach  number  and  time, 
respectively,  being  determined  either  theoretically  or  experimentally. 


The  traction  force  of  the  rocket  is  considered  constant  in 
the  fuel  combustion  stage  (active  stage)  and  it  is  assumed  that 
it  is  linearly  decreasing  in  time  to  zero  in  the  short  transition 
period  from  active  to  passive  stage . 
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3.  THE  DIFFERENTIAL  EQUATIONS  FOR  THE  MOTION  OF  THE 
ROCKET 


It  is  assumed  that  the  flight  distances  are  short  compared 
to  the  radius  of  the  Earth  so  that  the  acceleration  of  gravity  is 
considered  to  be  constant  in  magnitude  and  direction  (g  =  const.) 
and  that  the  rocket  motion  takes  place  on  a  linear  trajectory  tilted 
at  the  angle  ©  from  the  horizontal  direction,  situated  in  the 
vertical  plane  defined  by  the  direction  of  the  initial  velocity  of 
the  center  of  mass  of  the  rocket  and  the  vector  of  the  acceleration 
of  gravity  (fig.  1) . 

The  differential  equations  of  the  motion  of  the  rocket  in 
this  plane  are: 


T  cos  x  —  JD  —  mg  sin  6  =  mi, 

T  sin  a  -f  L  —  mg  cos  0  =  0, 
J/„  —  $(£  cos  x  D  sin  x)  »  J* . 


(l.a) 

(l.b) 

(l.c) 


In  these  equations  the  mass  of  the  rocket  m  has  a  linear 
variation  with  respect  to  time  in  the  active  stage,  being  constant 
in  the  passive  stage 


m 


m«(l  -Qt)  for  te  [0,0. 
m,  —  const.  for  t  >  i„ 


(2. a) 
(2.b) 


where  the  adimensional  flow  of  burnt  gases  is  constant  and  given 
by: 


Q  * 


(3) 


The  traction  force  T  is  considered  constant  in  the  active  stage, 
and  in  the  transition  period  it  decreases  continuously  in  time 


to  cancellation: 


T 


Tt  =  const,  for 
•  Ttt(t)  for 

0  for 


t  6  [0,  tt  ftr]« 

f  e  [f*  —  tfn 

t  >  t.. 


(4) 


A  linear  variation  of  the  traction  force  in  time  for  the 
transition  period  was  considered  in  the  numerical  example  at  the 
end  of  this  paper: 


T,r(t) 


t-JJ.  -  W 


(5) 


The  aerodynamic  forces  and  moments  are  given  by: 

L  =  ^P7‘5Ct ;  CL  =  «C2  4- 
Z>  =  _L  ?T*SCa, 

Mm  =  ypF’SiC.;  <7.  -  «C1  +  +  «<£, 

where  the  drag  coefficient  is  given  by: 


\Cd,  +  for  6  [0,  t,], 

l Co,  -f-  -h  C*/i  for  >  *■» 


(6. a) 
(6.b) 
(6  . c) 


(7. a) 
(7.b) 
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where  is  the  base  pressure  drag  coefficient  which  appears 

when  the  engine  ceases  to  function. 

In/v(6)  the  other  derivatives  of  the  aerodynamic  coefficients 
(like  Cl)  were  assumed  negligible.  Also,  considering  a  canard 
configuration  for  the  rocket,  with  the  command  surfaces  placed 
in  the  front,  the  effect  of  the  vortex  created  by  them  on  the 
wings  could  be  neglected. 

The  aerodynamic  coefficients  and  their  derivatives  are  known 
functions  with  respect  to  Mach  number.  Thus,  we  assume  their  values 
known  for  CD%(Mi),kJMi)yC0„(Xt),CUMi)y<rm(Mi),  CKJf,),  C&Jtt)  at 

different  Mach  numbers  Jfit  i  e  {1,  2,  3,. . **},  situated  in  the  range  of 
velocities  of  flight. 

Also  assumed  known  are  the  values  at  different  moments  t^  of 

the  distance  5  between  the  center  of  the  mass  of  the  rocket  C  and 

the  fixed  point  F  with  respect  to  which  the  aerodynamic  rolling 

moment  M  is  determined,  as  well  as  the  values  of  the  moment  of 
m 

inertia  J  of  the  rocket  with  respect  to  the  rolling  axis  that 
passes  through  the  center  of  mass,  that  is  £(*«)  and  J(t^)  are  known 
at  different  moments  t^  ,  i  e  {1,2,3,. . n,},  during  the  fuel  combustion 
stage  [o,  tal . 

In  order  to  determine  the  values  of  the  aerodynamic  coefficients 
and  of  their  derivatives  for  a  Mach  number  M  as  well  as  the  values 
of  the  moment  of  inertia  J  and  the  distance  5  which  defines  the  posi¬ 
tion  of  the  center  of  mass  at  a  time  t,  a  numerical  method  of  the  Sline 
type  has  been  used. 

Atmosphere  density  p  and  sound  velocity  in  the  atmosphere  are 
variable  with  respect  to  the  altitude  of  flight  h  of  the  rocket  and 
are  given  by 


fc  =  /t0  -r  j  sin  6,  (8) 
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in  which  h  is  the  altitude  at  take-off. 
o 

Substituting  in  the  system  of  differential  equations  (1)  the 
expression  of  the  mass  m  given  by  (2. a)  and  the  expressions  of  the 
aerodynamic  coefficients  given  by  (6)  and  (7. a)  and  those  of  the 
aerodynamic  forces  and  moments,  this  system  becomes: 

T cos  *  — i-p F*S(Co.  +  kD*s)  -  m0(  1  -  Qt)  g  sinQ  =  m0  (1  -  Qt)  x,  .  (9  .a) 

2 

T  sin  a  +  —  pF2S( aCl  +  )  —  mn(l  —  Qt)  g  cos  0  =  0,  ( 9  .  b ) 

^-PF*57(aOi  +  ?0i  +  qCt)  -  i  P^5[(«CI  +  fJCi)  cos  at  + 

2  2 

-f  (Ca,  +  Bin  a]  =  Ja.  (9  .  c) 

4.  INTEGRATION  OF  THE  DIFFERENTIAL  EQUATIONS  OF  MOTION  OF 
THE  ROCKET 


The  system  of  differential  equations  (9)  is  a  nonlinear  system 
whose  integration  is  possible  only  by  numerical  methods.  For  this 
purpose  we  obtain  the  expression  of  8  from  (9.c): 

P - - -  ;  — - «  -  otfoi - ^<71  cos  «)- 

cos.  y?FW  v  1  / 

I  (10) 

Substituting  this  expression  in  ( 9 . b )  and  using  the  notations 

T  _  m  q  cos0  T.  J 

T*  = - - - ;  Q*  =  - - ,J*=— - - #  (11) 

I.rw  i-o  V*S  —  pF*SI 


—  oV*S  - 

2  P  2 

the  following  equation  is  obtained: 


(T*  sin  *  +  *Cl  -  6*)  (^Ci  -  -j  Cl  cos  a^  4-  C\  - 
—  «^<7i--yC2cos«j-*  01  -  y  (C».  +  sm  *]  =  °’ 
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which,  together  with  (9. a)  is  a  nonlinear  system  of  differential 
equations  of  second  order. 


Then,  introducing  the  supplementary  variables  V  =  x  and 

•  ^ ^/i .  a 

q  =  a,  the  system  of  differential  equations^  of  second  order  becomes 
a  system  of  first  order 


?  =  ,[. 


i  =  r, 

T*  cos  a  —  (Cp,  + 


* 

—  r*  sin  a  —  aC2 


cos  6  —  sin  0  I  * 


i  _  1  r  G*  -  T*  sin  a  -  aCl  ( ^  l  \ 

9  Cl  -C,icos«j  + 

+  - -  Clcos  aj  a  -b  “7^1? - y-  (Cd,  +  fcjjad)  sin  *j 


(13. a) 
(13. b) 

(13. c) 

(13. d) 


This  nonlinear  system  of  differential  equations  of  first 
order  can  be,  in  principle,  numerically  integrated  using  an  adequate 
method  such  as  the  Runge-Kutta  method.  At  every  step  the  value  of 
the  attack  angle  8  can  be  determined,  by  using  (10). 


It  is  noted  though,  that  in  many  cases,  the  derivative  Cl 
of  the  drag  coefficient  with  respect  to  the  attack  angle  8  has  very 
small  values,  and  sometimes  people  tend  to  neglect  it,  which  leads 
to  great  errors  when  using  the  integration  method  above  mentioned 
(eg  is  in  the  equation  (13. d)  in  the  denominator)  and  even  making 

O 

it  impossible  to  use  the  respective  method  (in  the  case  that  c£ 
is  negligible) . 


For  this  reason,  a  procedure  which  can  be  applied  even  for 
very  small  values  of  the  derivatives  eg  and  even  for  eg  *  0  will 
be  shown  here. 


This  procedure  consists  of  using  at  the  same  time  the  Runge- 
Kutta  method  for  the  numerical  integration  of  the  system  of  dif¬ 
ferential  equations  (13. a)  and  (13. b)  and  a  method  with  finite 
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differences  in  order  to  integer  the  differential  equation  (12) 

which  is  equivalent  with  (13. c)  and  (13. d) .  Regressive  differences 

of  second  order  are  used  for  this  purpose  in  order  to  obtain  the 
derivatives  d  and  d*  which  appear  in  equation  (12): 


=  i>i*v  + 

*i  =  Ci*i  +  c., 


(14. a) 
(14. b) 


For  the  integration  step  i  =  1,  at  which  only  the  values  of 
a  marked  i  -  1  and  for  the  step  i  =  2  for  which  only  the  values 
marked  i  -1  and  i  -  "2  are  known,  the  values  of  b^ ,  b2,  c^ ,  C2, 
from  (14)  can  be  calculated  with  the  following  formula  * 


6,  =  -—?—;  bt  «  — i— 

'-3a,.,  -  At  •*,_!  + 

2  ,  ..  .  v 

(16 .a) 

2  At  2 -At 

2 

Cj  — y  Cj  “• 

(At)*  * 

- -(*i-i  +  At*  *1-1). 

(At)* 

(16 .b) 

It  is  noted  here  that  the  initial  conditions  oig ,  ao  "do  are 
known. 


From  (12)  the  next  relation  can  be  obtained  with  which  the  value 
ou  of  the  corresponding  integration  step  i  can  be  calculated: 

(T*  ain  *,  +  C2«,  —  0*)  -  y  Ci  cos  a,  j  + 

+  01  J^*(c1a<  +  c,)  —  ^(7i - ~Cl  cos  a,j  a,  — 

~  "^7  *«  +  &t)  4-  ~  (Co,  +  &!>■?)  sin  *=  0. 

*  The  equations  (16)  are  obtained  by  setting  ai~2'  ai- 3  as  functions 
of  ^i-l»’^i_]_  from  the  formula  with  regressive  differences  of 

second  and  first  order  which  yield  a^_^,*dj  .  respectively  by 
substituting  the  result  in  the  equations  (157. 
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In  order  to  determine  the  unknown  value  of  dj_  from  the  tran¬ 
scendental  equation  (17)  a  method  of  successive  approximations  is 
used.  For  this  purpose  the  equation  (17)  is  written  in  the  form 


*<  =  /(*<), 


(18. a) 


where 


f(t  )  =  ~  —  C08  g|  j  ^  ~  Caij)  ( a,  —  sin  a<) 

D  —  (A  —  B  cos  **  -r  Ca,J) 


(18. b) 


.4  =  T*Cl  -f  ^ 
Z 

(19. a) 

B  =  -7-  T*Cl 

(19. b) 

c  -  kDCl 

■ 

(19. c) 

D  =  CiCl  -Cl| 

'ci-j'c+^c ib,y 

(19. d) 

E  =  G*Cl  -C*| 

r  z  *  \ 

d»,). 

(19. e) 

F  =  G*Ci. 

X 

(19. f) 

It  can  be  shown  that  for  the  usual  range  of  incidences  (-  12° ,+ 
12°),  that  is  (-  0.2,  +  0.2  rad)  and  for  usual  values  of  aerodynamic, 
gemoetrical  and  inertial  values  of  the  rocket,  the  function  f(ou) 
defined  by  (18. b)  satisfies  the  conditions 

sup  ;/«,(«,).  =■  X,  X  <  1, 

j  g  ^ _ q  o  x  Q  (20) 


In  these  conditions,  the  application  f  is  a  contraction  and  thus, 
according  to  the  Picard  -  Banach  theorem,  the  fixed  point  of  the  function 
f(a^),  that  is  the  root  of  the  equation  (17)  exists  and  is  unique,  and 


it  can  be  determined  using  an  iteration  method  with  the  help  of  the 
series  generated  by  the  relation: 


E  —F  cosa,,-!  -r  (A  —  B  cos  C's£,,-i)(«.-.m-i  — *in  \ 

I)  —  (A  —  Hcos  —  C ) 

It  can  be  demonstrated  that,  no  matter  what  is  within  the 
interval  (-  0.2  +  0.2)  rad,  the  series  (21)  is  convergent  towards 
the  fixed  point  and  the  error,  which  may  occur  d(a^,  )  if 

instead  of  the  root  a.  the  value  a.  „  is  taken,  is  given  by  the 
relation 


d( *i,  *,.„)  <  - *,.#).  (22) 

1  A 

As  a  starting  number  ct^g  in  the  iteration  to  determine  the 
angle  ou,  the  value  of  the  incidence  at  the  previous  step  can  be 
used,  that  is: 


(23) 

It  is  mentioned  that  by  writing  equation  (17)  under  the  form 
(18),  the  iteration  procedure  used  to  determine  the  angle  becomes 
very  rapidly  convergent  in  all  the  analyzed  cases. 

With  respect  to  the  value  cu  so  determined  the  attack  angle  of 
the  command  surfaces  can  be  determined  at  every integration  step. 

(The  values  for  <5^  *  q^and’a^  which  occur  in  (10)  are  given  by  (14)). 

a)  ESTABLISHING  THE  INITIAL  CONDITIONS  OF  THE  MOTION 

In  order  to  start  the  numerical  integration  procedure,  the  initial 
conditions  of  the  motion  have  to  be  known,  that  is,  the  following 
kinematic  parameters  at  launching  time:  h  ,  xn,  Vn ,  <*q.,  qQ . 


Assuming  that  the  launching  takes  place  from  a  fixed  ramp  or 
from  a  mobile  ramp  (i.e,  a  plane)  whose  motion  is  known,  then  the 
initial  parameters  hg,  Xg,  '  qg  =  olq,  are  known  as  well  as  the 
angular  acceleration  Sg  =  *dg. 

With  regards  to  the  initial  value  a0  and  these  must  be 
determined  in  such  a  way  so  that  the  rocket  follows  the  linear 
trajectory  from  launching.  It  is  mentioned  that  in  the  present 
paper  it  is  assumed  that  the  rocket  moves  on  a  linear  trajectory 
straight  from  launching.  In  a  future  paper  the  case  when  the 
rocket  first  moves  on  a  curved  trajectory  and  then  starts  a  linear 
motion  will  be  considered. 

To  determine  aQ  the  initial  conditions  of  the  motion  are  sub¬ 
stituted  in  (12) .  A  transcendental  equation  in  aQ  is  obtained 
which  is  solved  using  the  same  method  of  successive  approximations 
as  in  the  case  of  equation  (17) .  For  this  purpose  the  equation  is 
put  under  the  form 


«•  =“/«(«o)» 


(24. a) 


in  which 


^  )  =  -ge  ~  Fq  co*  «o  -4-  U,  —  B0  cos  gp  C0i$)  ( ap  —  sin  x„) 
A>  +  (.i#  —  S0  cos  *o  —  C0xi) 


(24. b) 


and  where: 
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It  can  be  demonstrated  that  in  this  case  also  fg(aQ)  is  a 
contraction  so  that  the  Picard  -  Banach  theorem  can  be  used  for  the 
convergent  series  of  successive  approximations: 


E9  —  F0  cob  gp.-i +(A0—B0  cos  ao..-I-f-<70ao,.-i)(a0  .-1— 8in  . 

D0  4-  (A%  —  B0  cos  4-  Ct  «o..-i)  ( 26  * 


The  starting  value  for  the  iteration  to  determine  aQ  is  obtained 
for  n  =  1  and  aQ.0  a  0  in  (26)  : 


*0.1  = 


(27) 


The  iteration  process,  which  is  rapidly  convergent,  stops  at  the 
number  n  which  yields  aQ  with  the  desired  accuracy  (aQ  ^  aQ  n)  . 


Having  aQ ,  the  initial  attack  angle  can  be  calculated  with 
(10),  in  which  the  initial  values  aQ,  dQ  =  q0,'dt0  =  eQ  are  sub¬ 
stituted,  as  well  as  the  initial  values  of  the  aerodynamic  coeffic- 

.  ... 

rents  of  other  parameters  variable  in  time,  such  as  m0,  £q»  J0'Aetc* 


b)  NUMERICAL  INTEGRATION  OP  THE  SYSTEM  OF  DIFFERENTIAL  EQUAT¬ 
IONS  OF  MOTION 

Knowing  the  initial  conditions  for  the  motion:  hg ,  Xq,  Vq ,  a0 ,  q0 
the  system  of  differential  equations  (13. a)  and  (13. b)  can  be  integr¬ 
ated  by  one  of  the  known  numerical  methods,  i.e  the  Runge  -  Lutta  - 
Gill  method  of  the  fourth  order,  which  has  the  desired  accuracy  using 
a  halving  method  for  the  integration  step  for  which  the  calculation 
error  is  reduced  to  the  imposed  value. 

The  angle  a,  of  (13,b)  is  determined  at  every  step  using  the 
iteration  procedure  above  mentioned. 


During  the  numerical  integration,  the  aerodynamic  coefficients 
and  their  derivatives  are  considered  variable  with  Mach  number  being 
given  by  numerical  values.  Their  values  corresponding  to  the  Mach 
number  of  the  flight  at  a  given  moment  are  calculated  using  a  Spline 
type  of  numerical  interpolation  method.  A  similar  method  can  be 
used  to  determine  at  every  integration  step  the  moment  of  inertia 
J  of  the  rocket  and  the  coordinates  £  of  the  center  of  mass  which 
also  are  variable  with  respect  to  time,  being  previously  calculated 
and  given  by  numerical  values.  They  can  be  calculated  also  during 
the  numerical  integration  at  every  step,  if  the  relation  between  the 
fuel  consumption  and  its  distribution  in  the  rocket  body  is  known. 

To  calculate  the  thermodynamic  parameters  of  the  atmosphere 
(density  and  sound  velocity)  which  are  variable  with  the  altitude 
h,  a  method  of  numerical  interpolation  can  be  used  starting  from  their 
values  in  standard  conditions  or  using  analytical  calculus  relations 
as  is  the  case  in  the  numerical  example -presented  at  the  end  of  this 
paper . 

At  every  step  of  numerical  integration  which  corresponds  to  a 
certain  value  of  the  time  t,  by  the  algorithm  presented  in  part  4 
the  values  of  x,  V,  a,  q,  e  are  determined,  with  which  from  (ID) 
the  value  3  of  the  attack  angle  is  determined. 

£ 

The  succession  in  time  of  the  values^ of  the  attack  angle  repres¬ 
ents  the  required  command  law  for  the  rocket  to  move  on  a  linear  tra¬ 
jectory. 

5.  THE  MOTION  OF  THE  ROCKET  IN  THE  PASSIVE  STAGE 

(t  >  t.h 

The  rocket  has  to  continue  the  linear  motion^ after  engine  shut-off 
The  differential  equations  of  the  motion  for  this  stage  are  obtained 
from  the  corresponding  equations  in  the  active  stage  in  which: 

T  =  0,  J  SB  Jf  am  const., 
m  “  *'i  (1  QU)  *»»/*■  const.,  5*5/=*  const. 
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and  the  bottom  drag  coefficient  is  added  to  the  drag  coeffic¬ 

ient  at  zero  lift 

Do 

For  the  transition  period  it  can  be  assumed  that  for  the  time 
ttr  =  (0-2  -  0.4)  seconds  the  mass  of  the  rocket  decreases  linearly 
in  time  as  in  the  active  stage,  and  the  traction  force  for  the  trans¬ 
ition  period  Ttr  also  decreases  linearly  from  a  value  Tq  to  zero 
according  to  (5) .  In  order  to  find  the  differential  equations  of 
the  motion  for  this  stage,  we  substitute  T  with  Ttr  in  the  equations 
of  the  motion  for  the  active  stage.  The  same  numerical  methods  will 
be  used  as  for  the  active  stage  for  integration  of  the  differential 
equations  of  the  motion  of  the  rocket  corresponding  to  the  last  two 
stages  (transition  and  passive) . 

6.  NUMERICAL  EXAMPLE  OF  CALCULATION 


In  order  to  illustrate  the  above  shown  algorithm,  the  following 
numerical  example  is  presented  in  which  these  data  have  been  used: 

-  geometrical  and  inertial  characteristics  of  the  rocket: 


5  =  1  *»,  =«  600  kg, 

l  —  0,7  m,  m,  *  400  kg. 


5f  =  1»3  m, 


</„  =*  124, ."3  kg  m*, 


5c,  =  1,86667  m,  Jf  =  100  kg  m*  ; 

5c,  =  1,8  m, 

-  the  characteristics  of  the  motor  ( three  successive  values 
have  been  considered  for  the  traction  force) : 

To  =  25  000  X,  t,  =  3  s, 

To  =  50O00X,  /„  =  0,4  s; 

T0  =  100  000  X, 


J 
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-  the  initial  conditions  of  the  motion: 


h0  =  1  000  m, 
x0  =  0, 


1'0  =  200  m/s, 

2o  =  «o  =  0 ; 


-  the  slope  of  the  linear  trajectory 

0  =  15°; 

-  the  required  time  for  the  rocket  to  maintain  a  given  linear 
trajectory 


tf  =  25  s. 

The  position  £c  of  the  center  of  the  mass  of  the  rocket  and  the 
moment  of  inertia  J  with  respect  to  the  -rolling  axis  have  been  cal¬ 
culated  as  numerical  functions  of  time  considering  a  certain  fuel 
distribution  inside  the  rocket.  Their  variation  in  time  is  shown  in 
fig.  2  and  has  been  introduced  in  calculations  as  numerical  form. 

The  variation  of  the  aerodynamic  coefficients  and  of  their 
derivatives  with  respect  to  the  Mach  number  is  indicated,  for  the 
example  given,  in  fig.  3,  4,  5.  Their  values  have  been  introduced 
in  calculations  in  numerical  form.  kD  Q  (1/rad)  and  CDfd  *  0 
were  considered  in  this  example.  Air  density,  variable  with  the 
altitude  h,  was  calculated  by: 

F  —  Pr'e~ip*i  (2  8) 

.  3 

where  pp  is  air  density  at  sea  level  (pp  =  1.229  kg/m  0) ,  and  kp  is 

a  constant  for  which  different  values  were  taken  in  1,000  m  ranges, 

values  calculated  with  the  data  derived  from  standard  atmospheric 

values. 
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In  order  to  calculate  the  sound  velocity  in  the  atmosphere 
the  following  expressions  have  been  used: 


—  Op  —  1ish 

for 

h  <  10  •  668  m, 

—  296,2656  m/s 

for 

h  >  10  •  668  m, 

(29) 


where  ap  is  sound  velocity  at  sea  level  {ap  =  340.1568  m/s),  kg  is  a 
constant  (kg  =  0.004114286  s-1)  and  h  is  measured  from  sea  level. 


The  results  are  shown  in  the  diagrams  of  fig.  6  -  10  in  which 
the  distance  x  that  the  rocket  cruised,  the  velocity  of  the  rocket 
V,  the  Mach  number  of  the  flight  M,  the  angle  of  incidence  a  and 
the  attack  angle  B  of  the  command  surfaces  were  represented  with 
respect  to  time  for  the  three  values  considered  for  the  traction  force. 


Fig.  2 
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Fig.  10 


The  variation  in  time  of  the  attack  angle  3  of  the  command 
surfaces,  shown  in  fig.  10,  is  the  command  law  for  the  rocket  to 
follow  its  linear  trajectory  in  the  numerical  example  presented  here. 
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